Hydromagnetic Energy Conversion and Prompt Collimation 
in Mildly Advective, Kerr Black Holes* 



Rafael A. Araya-Gochez 

Laboratorio de Investigaciones Astrofisicas, Escuela de Fi'sica 
Universidad de Costa Rica, San Jose, Costa Rica. 



o 
o 
o 

(N 

c 

c3 



> 

-xf : 

m ■ 

(N ■ 

On : 
ON ■ 

Of 

6 



Abstract 

Recent evidence of the phenomenal energetics involved in 7- 
ray burst events strongly suggests that the progenitor must 
efficiently convert gravitational binding energy into a moder- 
ately collimated outflow, possibly in the form of a Poynting 
jet. We show that an MHD-instability driven dynamo (IDD) 
operating in a hot accretion disk is capable of generating 
energetically adequate magnetic flux deposition rates above 
and below a mildly advective accretion disk structure. The 
dynamo is driven by the magnetorotational instability (MRI) 
of a toroidal field in a shear How and is limited by the buoy- 
ancy of 'horizontal' flux and by reconnection in the turbulent 
medium. In the comoving frame of a semi-thick, slim disk 
setting, the predominant field components reside in surfaces 
perpendicular to the local meridian and the flux is deposited 
in funnels that are relatively free of baryons. The efficiency 
of magnetic energy deposition is estimated to be comparable 
to the neutrino losses but the strong effective shear induced 
by the metric on the MRI favors pumping magnetic field 
energy at low wavenumbers, i.e., field generation at large co- 
herence lengthscales. This, in turn, suggests that an MHD 
collimation mechanism may deem this process a more viable 
alternative to neutrino-burst-driven models of 7-ray bursts. 



1 Introduction 

The combined redshift and fluence measurements of at least 
five 7-ray burst sources: GRB970508 @ z = .835 (Metzger 
et al. 1997, Kouveliotou et al. 1997), GRB971214 @ z = 3.4 
(Kulkarni et al. 1998, Kippen et al. 1997), GRB980718 @z = 
.966 (Djorgovski et al. 1998, Kippen et al. 1998), GRB990123 
@ z = 1.6 (Kelson et al. 1999, Kippen 1999a, Conners 
et al. 1999), and GRB990510 @ z = 1.619 (Vreeswijk et 
al. 1999, Kippen et al. 1999b), plus very large photon energy 
detections in certain bursts (e.g. Sommer et al. 1994, Hurley 
et al. 1994, Dingus 1995) and tight size constraints derived 
from the rapid risctimes of burst triggers (Walker, Schaefer 
& Fenimore 1998) strongly suggest that the release of energy 
is highly focused by the central engine that propels a 7-ray 
burst. 

In spite of the very large 7-ray energy requirements, E 7 ~ 
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( x I.O+54 erg (£ = <5fi 7 /47r), the efficiency of energy deposi- 
tion into electromagnetic channels is likely to be very poor, 
Ej/Etot = e s[.01, .001], if the burst is driven by a neu- 
trino burst in analogy to the processes thought to give birth 
to supernovae (MacFadyen & Woosley 1998, 1999) or if the 
burst involves major energy losses to gravitational radiation 
such as might be the case in compact object merger scenarios 
(Rasio and Shapiro 1992, 1994, Davies et al. 1994, Janka k, 
Ruffert 1996, 1998; Ruffert & Janka 1998). 

Thus, although the monumental observational progress 
of late at first prompted some workers to suggest that 
"the physical mechanisms behind 7-ray bursts (were) within 
reach" (Metzger et al. 1997), the unexpectedly large increase 
in the energy and/or collimation requirement of any viable 7- 
ray burst progenitor has deemed such mechanisms at present 
largely undetermined. Indeed these measurements pose a se- 
rious energy budget problem for arguably all gravitational 
collapse powered theoretical models of 7-ray bursts if the 
energy release is not moderately collimated £s 1.0_ 3 (see, 
however, Meszaros, Rees & Wijers 1999). 

On the other hand, the angular spreading problem (Feni- 
more, Epstein & Ho 1996, Fenimore, Mandras & Nayakshin 
1996) means that it may be quite implausible for complex 
millisecond 7-ray variability to be attributed to the interac- 
tion of a single fireball with an external medium. The some- 
what favored scenario is one in which the variability down 
to possibly sub-ms timescales portrays flux fluctuations at 
the emission radius R e (Sari & Piran 1997, Kobashi, Piran 
& Sari 1998) or from flaring within the energy source itself. 
Fenimore, Ramirez- Ruiz & Wu (1999) assert that the latter 
must be the case for the brightest burst yet: GRB990123. It 
is not clear, however, how such radiation could avoid ther- 
malization unless the energy release mechanism yields in- 
finitesimally thin shells (M. Baring, Priv. Comm.). 

An attractive solution to the energy budget/collimation 
problem starts with a directional Poynting-flux dominated 
outflow (Thompson 1994, Meszaros & Rees 1997) under the 
premise that such a flow may carry very little baryonic con- 
tamination if deposited along a centrifugally (or gravitation- 
ally) evacuated funnel such as the angular momentum axis 
of a black hole-accretion disk system. 

The notion that magnetic fields may play an important role 
in tapping the energy available in accretion disks and in its 
subsequent re-processing into a high Lorentz factor outflow 
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has been surmised by several authors. Attempts have been 
made at drawing an analog to the standard phenomenol- 
ogy of pulsar electromagnetic emission (e.g, the magnetized 
differential rotor Katz 1997) and at invoking the Blandford- 
Znajek mechanism (BZm) (Paczyhski 1998) and BZm-like 
processes at the innermost regions of accretion disks (Ghosh 
& Abramowicz 1997, Livio, Ogilvie & Pringle 1998). 

These estimates of Poynting flux luminosities assume a pri- 
ori the existence of a nearly equipartition, large-scale (0[r\), 
external magnetic field especially with regards to jet forma- 
tion processes (Blandford & Payne 1982, Tsinganos & Sauty 
1994, Scheuer 1994, Blandford 1994) even if the magnetic 
collimation mechanism does not explicitly involve significant 
baryon entrainment (Lyndcn-Bcll 1997, Lynden-Bell & Boily 
1994, Appl & Camenzind 1993). Yet, such an assumption is 
hard to justify on theoretical grounds. On the one hand, 
self-generated magnetic fields in thin Keplerian disks are un- 
likely to attain coherence lengthscales far in excess of the 
disk's pressure scale height, TL, and inverse-cascade estimates 
in thin disks (Tout and Pringle 1995) indicate that the large 
scale component is much weaker 0{TL/r). On the other hand, 
advection of a large scale fields from the outer disk or stellar 
mantle is unlikely to permit a substantial steady inflow of 
material (Ghosh & Abramowicz 1997) and thus preclude the 
large accretion rate required by the overall energetics. 

In addition, jet-launching by large scale magnetic fields is 
known to possess a large degree of baryon entrainment (see, 
e.g. Pringle 1993). This 'problem' (from the 7-ray burst 
standpoint) persists even when the collimation mechanism 
does not invoke the presence of a baryonic wind (e.g. Lynden- 
Bell & Boily 1994, Lynden-Bell 1997) because a large scale 
field connects baryon-free with baryon-loaded regions of the 
system. Thus, just how the Poynting jet might form is an 
outstanding theoretical problem posed to remain a subject 
of controversy and future research for some time. 

On the other hand, to our knowledge formal motivation for 
an external field with the desired strength has yet to be inves- 
tigated in this setting (with the possible exception of Thomp- 
son (1994) who invokes magnetar-like processes, i.e. strong 
neutrino induced convective instabilities at near nuclear den- 
sities). Simple 'winding-up' of radial field to produce a 
toroidal component was invoked by Narayan, Paczyhski & 
Piran (1992) and by Katz (1997) but this mechanism does 
not replenish the radial field thus failing to operate as a dy- 
namo. Parker (P-type) and magnetorotational instabilities 
(A4-type or MRI's for short) are also often invoked with- 
out further elaboration on how this processes might operate 
together under the conditions implied by 7-ray burst mod- 
els (see, however, Tout & Pringle 1992, hereafter TP92, for 
a general, quantitative account of an instability driven dy- 
namo). 

We envisioned the accretion disk setting following the stan- 
dard hyper-accreting black hole model of Popham, Woosley 
& Fryer (1999, hereafter PWF): M hh = 3 M ; a SS = 
0.1, and M = 0.1 M Q sec -1 . These authors find that the 
onset of photodisintcgration and of neutrino cooling at radii 



r £ 70 [GM/c 2 ] yield a mildly advective accretion disk struc- 
ture]] (Tt@/r s.4) with semi-Keplerian rotation rate. The 
pressure at r ~ 70 is dominated by the nucleon gas compo- 
nent and inside this radius, it becomes dominated by radia- 
tion at r s 40. The temperature rises above (I7) pair creation 
threshold at r ~ 56. For a viscosity parameter ass = -1, 
a fluid element's total number of windings is very limited, 
^Vwind ~ 10, but this increases almost tenfold by reducing 
&SS by the same factor. The total neutrino luminosity, about 
2% of Mc 2 , and annihilation rate efficiency, « 12%, for this 
model are probably insufficient to explain moderately bright 
7-ray bursts but the authors remark that magnetohydrody- 
namic processes might be more efficient at producing a high 
energy (Poynting) jet. 

We elaborate on this suggestion using an instability-driven 
accretion disk dynamo (IDD) mechanism. Similar versions 
of this model (Tout & Pringle 1992, hereafter TP92) involv- 
ing V-type and axisymmetric A4-type (i.e. Balbus-Hawley) 
instabilities have meet with limited success in explaining 
dwarf novae quiescence and eruption phenomenology (see Ar- 
mitage, Livio & Pringle 1996 and references therein). How- 
ever, invoking Parker's undulate instability to promote the 
growth (loss) of vertical (horizontal) field is questionable in 
turbulent disks where the turbulence is fed by MPJ insta- 
bilities (Vishniac & Diamond 1992), unless the growth and 
length scales of P-m odcs could compare favorably to the 
MPJ scales (see §^8). It is at present unclear what role, 
if any, V modes may play in the present context where the 
(presumably) thermal pairs are collisionally well coupled to 
the nucleon fluid. Following this argument, 'P-modes are as- 
sumed to be suppressed thus precluding the meridional field 
from becoming dynamically significant. By contrast, MRI's 
will develop at the largest lengthscales permitted by the lo- 
cal field topology in spite of the turbulent component to the 
flow/field. The field lines should be, in fact, only mildly 
stochastic on the largest MPJ lengthscale. 

The dynamics of a dynamo process occurring inside an 
accretion disk are likely to be fairly sensitive to the equilib- 
rium field structure and topology. Because MRI's are surely 
present in any highly conductive astrophysical disk, no well 
defined global field topology is expected (inside the disk). In 
addition, the field structure could be either diffusive or in- 
termittent, i.e. concentrated in semi-empty flux ropes. In 
a recent paper (Araya-Gochez 2000a), we have looked into 
the case of an intermittent field in a static metric. Here we 
consider both cases in a more general stationary metric. 

The field components correspond to azimuthally averaged 
quantities in the comoving frame, i.e. in the local standard 
of rest frame. These are consistent with radial field gen- 
eration from the MRI evolution of a toroidal magnetic field 
(i.e. non-axisymmetric Al-type instability) which is, in turn, 
generated by the sheared radial field; with both components 
limited by the buoyancy of nearly 'horizontal' flux and by 



1 Here the scale height, "He = rSfd corresponds to an average on 
spherical surfaces consistent with a not so thin, partially advective disk 
(see, e.g., Abramowicz, Lanza & Percival 1997). 
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fast turbulent reconnection of field lines. This model is only 
weakly sensitive to the details of the turbulent state but it 
does depend explicitly on the magnetic Mach number of the 
turbulence. 

Although MHD disk simulations with vertical stratifica- 
tion generally fail to support the buoyancy of magnetic fields 
(Brandenburg, Nordlund, Stein & Torkelsson 1995; Stone, 
Hawley, Gammie & Balbus 1996), these results prove to be 
very sensitive to the choice of vertical boundary conditions 
and do not account for the effect of radiation/pairs (see §B 



Lastly, we note that the field estimates are not explicitly 
sensitive to the precise value of the viscosity parameter, ass, 
but depend rather on local pressure ratios, on the relativitic 
generalization of the shear parameter (A{^' rt = Oort's A con- 
stant), and on the linear timescale for development of the 
non-axisymmetric MRI (which is found to differ from Aq^ 4 
by a general relativistic correction factor, §A). Angular mo- 
mentum may be transported by the ordered and/or the tur- 
bulent components to the field but the relative contributions 
to the viscosity are not clear at present. 

A review of the slim disk setting in the Kerr geometry and 
associated disk thermodynamics is given in §|2| The instabil- 
ity driven dynamo in a mildly advective disk is laid out in 
§|j| We then critique and review the assumptions, find steady 
state solutions to the dynamo equations and use these to esti- 
mate efficiency of the hydromagnetic energy deposition from 
field buoyancy in §|[ 



2 Accretion Disk Setting 

2.1 Slim disk in Kerr geometry 

The Boyer-Lindquist generalization of the Schwarzschild co- 
ordinates (i, r, 9, ip) affords the most popular form of the 
metric for a rotating Kerr black hole. In normalized geomet- 
rical units (c = G = Mbh = 1), the angular momentum of 
the hole is a = J [l/(GM hh /c 2 )}. 

The metric functions of the radial BLF coordinate are 
written as relativistic corrections (e.g. Novikov and Thome 
1973)0: 

A = l + a 2 /r 2 + 2a 2 /r 3 , and V = 1 - 2/r + a 2 /r 2 , (1) 

so that, in the equatorial plane, the squared line-element may 
be written as: 

ds 2 = -^dt 2 + r 2 A{d4> - uidt) 2 + ^dr 2 , (2) 

with uj = 2a/ Ar 3 the rate of frame dragging by the hole. 

At 9 = it/2, the horizon lies at r + = 1 + \/l — a 2 (V\ r± = 
0) and no observers with time-like worldlines may keep ip 
constant inside the static limit at rs = 2. 

The Killing vector fields that belong to the azimuthal and 
time symmetries of the metric arc deliberately simple in the 

2 A useful relation connecting A and 23 is: (1 — 2/r)^t+4a 2 /r 4 = 23. 



BLF: * = (0,0,0,1) and T = (1,0,0,0). These may be 
used to show that the time and azimuthal components of the 
four- velocity, (U t ,U r ,Ug,U v ), are conserved along geodesies. 
We denote the latter according to standard notation: Ut = 



Pi 



-8, and U w 



The unit time-like vector is constructed from the killing 
vectors: 1( cx T + CijM/. This expression defines "forwards in 
time" everywhere in the flow. The rate of frame dragging and 
the radius of gyration for circular orbits are quantities de- 
rived from inner products of Killing vector fields (Abramow- 
icz, Chen, Granath & Lasota 1997): w = -(Y*)/(**) and 
zu = (\E"J>)/|Y respectively. 

A natural frame to treat the local aspects of disk flow is 
the locally non-rotating (LNR) frame which is dragged by 
the hole at a rate u> with respect to distant observers and 
where l t - = (1, 0, 0, 0). Bardeen, Press and Teukolsky (1972) 
pointed out the advantages of calculating physical processes 
there instead of in the Boyer-Lindquist frame which is fixed 
with respect of distant stars. However, in the spirit of the 
shearing sheet approximation used in the next section, the 
frame we choose to evaluate the ingredients necessary for dy- 
namo action is the local rest frame (LRF) of the fluid where, 
in general, l t - =^ (1, 0, 0, 0) (in our notation t denotes the time 
component measured in the LNR frame while t denotes the 
same measured in the LRF). 

The flow moves with respect to the LNR frame with az- 
imuthal three speed v v = x v /x l (= (3 V ). This is related to 
the angular velocity, Q, as measured by distant observers by 



zu x (fi — lu) = r 



A 



V 



x 1 



r7V A 



(3) 



From the frame that corotates with the fluid, the measured 
radial velocity is defined as (3 r . This choice of frame to mea- 
sure the radial velocity dictates that j3 r < 1 everywhere in 
the flow (reaching the upper limit at r+) and that j3 r w c s 
at the sonic surface (Abramowicz et al. 1997). 

The orthonormal thetrad of basis LRF vectors, used to 
transform BLF <^=> LRF, and explicit forms for the four ve- 
locity as function of [3 r and t in the BLF are given in Gammie 
and Popham (1998, hereafter GP98). For completeness we 
exhibit the contravariant components (i.e. the momenta) 
which are easily derivable from the above definitions and 
from V/p u = —7: 



(U t , U r , U e , U v ) = 



(4) 



where 7 = 7 v 7 r , with 7^ and 7 r as measured in the frames 
indicated above. 

By symmetry and for simplicity U e is assumed to vanish at 
the equatorial plane and periodic azimuthal boundary con- 
ditions in the "dragged" frame allow for non-axisymmetric 
unstable MHD modes to develop while permitting axisym- 
metry of the mean flow. Although the equations governing 
the flow are greatly simplified by the symmetry conditions, a 
potential caveat of this assumption is that low wavenumber 
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non-axisymmetric unstable azimuthal modes, such as undu- 
late "P-typc modes, may violate these conditions if the mean- 
fields become kincmatically important. 

In accordance with the standard (one-dimensional) slim- 
disk approach (Abramowicz et at 1989, Narayan & Yi 
1995ab, GP98), the 'vertical' average is performed over the 
meridional and azimuthal coordinates 



(5) 



Tin is the ^ _ 



where: / is any dynamical variable, O-h = 
meridional angular pressure scale height of the disk (formally 
defined below) and where 



2,2 2 

a cos i 



.9 = |Det(g Mt ,)| = (r 
is the determinant of the full metric 



r 2 sin 2 i 



With T» v = pg^ + qu»u v 



Srisc' * ne y i scous stress-energy 



tensor, four vertically averaged equations for I, (3 r , T, p 
and for the angular momentum eigenvalue, j, may be derived 
from particle number, (pu^)-^ — 0, and energy-momentum, 
T^ v . v = 0, conservations (GP98, PWF, Abramowicz et 
al. 1997). In §A, we will use the radial momentum con- 
servation equation 



dfi r 
dr 



r 2 V 



1 - 



n 



1 - 



g dr ' 



(6) 



where p is the total pressure (including the magnetic contri- 
bution, e.g. §2.2), g = p+u+p is the (unnormalized) specific 
enthalpy and Q± — ±(r 3 / 2 ±et) -1 are the prograde and retro- 
grade Keplerian rotation frequencies of circular, planar orbits 
as measured in the BLF. 

PWF find that cooling from neutrino losses and from pho- 
todisintegration of He nuclei becomes important as the tem- 
perature and density rise. The cooling is calculated explic- 
itly (i.e. by solving the global ADAF). The disk is found to 
be optically thin to neutrinos for Msl.O M Q sec -1 , so that 
neutrino producing processes, pair annhilation and pair cap- 
ture on nuclei, act to cool the medium at moderate accretion 
rates. 

Where the mass fraction of nucleons X nuc < 1 , photodis- 
integration quickly cools off the disk. At moderate accretion 
rates, MzO.l M Q sec -1 , photodisintegration is complete at 
radii r s 70 0. The nucleon component may then be assumed 
to behave like an ideal gas; decoupled from the radiation field 
except to maintain an equilibrium temperature with the ther- 
mal bath (no two-temperature model is expected as in the 
optically thin case). 

2.2 Fluid component coupling 
and effective sound speed 

Ignoring electron degeneracy pressure, the equation state for 
the medium has two separate contributions p = pT + a'T , 

3 The nucleonalization radius should move farther out with the spin 
parameter, a, of the hole as it tends to increment the density at constant 
temperature (see next footnote). Evidently, higher accretion rates move 
this radius outwards as well. 



where a' = y§ a rad accounts for radiation and pairs. The 
effect of degeneracy in the equation of state is found (by 
PWF) to be minimal for M < 0.1 M Q sec -1 . Thus, unless 
the cooling is significantly enhanced; e.g. if the dynamo con- 
stitutes a significant sink of energy from the shear flow, one 
may cautiously^ neglect electron degeneracy in the disk. 

With the above approximation, the nucleon-gas+radiation 
medium is expected to posses an effective adiabatic index 
ranging from that of a pure relativistic radiation-e pair 
gas to that of a non-relativistic, ideal, heavy, nucleon gas 
r = (3[i np ] lnp) s £ [4/3, 5/3]. The pressure from an isotrop- 
ically tangled (i.e. turbulent) magnetic field is included by 
noting that such field couples well to the (heavy) nucleon-gas 
component through flux freezing, even though the turbulent 
field component by itself would be expected to behave like 
a relativistic gas (e.g. Narayan & Yi 1995b, Quataert & 
Narayan 1999). The coherent and turbulent component con- 



tributions to the pressure are p^ 



cohc 



turb 

Pb 



Pb = b /24-7T respectively, with B 



local 



B 2 /8tt, and 
B + b. 



Defining j3 = p Ta ,d/ (Pgas + Pb), and assuming a relativistic 
thermal pair distribution function, the relativistic sound 
speed is simply (Chandrasekhar 1939, Mihalas & Mihalas 
1984) 



with r 




(7) 



Conversely, were there a significant fraction of pairs populat- 
ing suprathermal energies in a less dense environment, the 
former would couple to the mean-field component through 
Parker's (1965) suprathermal MHD mo de t hus yielding two 
natural sound speeds for the medium (§3.6 



2.3 Vertical structure 

In the absence of detailed knowledge on the 'vertical' depo- 
sition of shear energy, little can be said about the thermal 
gradient along lg. The entropy gradient is thus assumed to 
be convectively stable. Strong stratification, as it occurs in 
neutron stars where the field has to reach super-equipartition 
values to break free buoyantly (Thompson 1994), is unlikely 
to occur in a disk where self-gravity is negligible. 

Narayan and Yi (1995a) have pointed out that, at least in 
the self-similar regime, ADAF's are well described by "ver- 
tically" integrated versions of the fundamental equations as 
long as the boundary surfaces are spherical in flat geometry. 
The averaging prescription of Eq || (note a typo in GP98), 
generalizes this procedure to curved three-space. Because of 
curvature effects, He = rQu J y/gee dd, i.e. the pressure 
scale-height as calculated in the LRF does not correspond to 



4 PWF find that when the rotation parameter of the hole is large, 
a ~ .95, this results in substantial enhancement of the density at the 
innermost region of the disk (rSlO). An accompanying increase in the 
efficiency of neutrino cooling then acts as a thermostat precluding the 
temperature from rising. Whether this situation makes the degeneracy 
pressure contribution to the equation of state non-negligible is unclear 
though. 
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a locally integrated meridional line clement. tig is merely a 
convenient definition of "height" using the BLF V coordi- 
nate. This definition breaks down when the disk is very thick, 
Qu > 1, but for 8-HS.4 the correction is sma ll (w 2.5% for 
a = .95). The instability lengthscales in §3.4 should corre- 
spond to locally integrated line elements but, ignoring this 
small correction, the pressure scale height is used as the stan- 
dard normalization. 

Abramowicz, Lanza and Percival (1997, hereafter ALP) 
have derived an equation for the pressure scale height which 
contains only thermodynamic potentials, geodesic flow in- 
variants, and the meridional 4- velocity component and its 
radial gradient. In the present situation, buoyant instabil- 
ities may render the velocity terms non-negligible, but one 
may still expect the fluid to be subsonic or transonic far from 
the horizon; effectively "freezing out" close to the horizon 
(i.e. flowing along 9 ~ const.). In this case, ALP give 



rin 



po i_ 



1< 

r m 



(8) 



where C — £ 2 — a 2 (£ — 1), c s is the effective soundspeed at 
the midplane of the disk and Q z is an effective frequency of 
meridional oscillatory motion (GP98). 



3 An Instability-driven Dynamo 

In spite of being the subject of great theoretical effort since 
its inception (Sakura & Sunyaev 1973), clear resolution on 
the question of the origin of anomalous viscosity in accre- 
tion disks (Pringle 1981) has only recently occurred. Hawley, 
Balbus and Winters (1999) conclude that non-magnetic hy- 
drodynamical Keplerian accretion is stable (see also Godon 
and Livio 1999) and that shear-fed MHD turbulence in a 
fully developed state is capable of explaining anomalous an- 
gular momentum transport in accretion disks. (Balbus and 
Hawley 1991 (BH91), 1992a (BH92a), 1992b, 1996, 1998; 
Hawley and Balbus 1991, 1996; Hawley, Gammie and Bal- 
bus 1995; Brandenburg, Nordlund, Stein & Torkelsson 1995, 
1996; Stone, Hawley, Gammie & Balbus 1996; Brandenburg 
1998). 

MRI induced magnetohydrodynamic turbulence may thus 
be reasonably proclaimed to be the source of anomalous an- 
gular momentum transport in accretion disk systems and 
the question of whether large scale fields are present must be 
subjected to their co-existence with the turbulent flow. 

Because the IDD is shear-fed, where the disk is advective 
and gas pressure dominated the low shear rate of an ADAF 
composed of non-relativistic particles gas starves the dynamo 
out (recall that when r e ff ec tive —> 5/3, self-similar solutions 
have nearly spherical accretion O — ^ 0, Narayan & Yi 1994, 
Narayan, Mahadevan & Quataert 1998). We thus consider 
the case where a large-scale field in the outer disk begins to 
emerge where its influence is dynamically unimportant (e.g., 
at r w 70, once photodisintegration is complete). 



3.1 On the effect of a radial flow 

Although a predominantly toroidal field may be expected a 
priori (at least on surfaces perpendicular to lg), the radial 
flow is non-negligible in the innermost disk. For instance, in 
the standard model of PWF, the ratio SJ v/ y = ("f^ 1 Pr) / Pip 
which compares the radial and azimuthal three- velocity com- 
ponents as measured in the LNR frame, is small but signifi- 
cant in the range r ~ [70, 4.5]. At r ~ 14 (where 9-h ~ -4 
is maximum), the sub-Keplerian (3 V ~ .26 and 5R ~ .12. At 



4.5, & 



.62 and 3? ~ .32 while at 



3.9, fa 



.61 



and -ft ~ .43 and quickly rising. 

We envisioned that the MHD condition of high conduc- 
tivity in a strongly sheared semi-Keplerian flow produces a 
mean-field on 9 = constant surfaces that tracks the mean- 
flow on average (equivalently, one may assume vanishing an- 
gular velocity gradients along mean field lines). Thus, since 
the presence of a small radial field component does not alter 
the Balbus-Hawley instability (BH91), one may cautiously 
ignore the effect of a small, passively advected, B r compo- 
nent also on the (linear) development of non-axisymmetric 
.M-type instabilities as calculated form a purely toroidal field 
( §|3.7j ). Note, however, that these two MRI's originate from 
different branches of the MHD wave dispersion relation (the 
slow and Alfven branch). This assumption could require re- 
evaluation when B r zB v . 

3.2 On the use of spherical coordinates 

To date, most magnetohydrodynamical stability analyses of 
accretion disk systems have been performed in the context of 
the thin disk formulation where the use of cylindrical polar 
coordinates is customary. In a semi-thick, slim disk setting, 
the polar z coordinat e is preferably replaced by the merid- 
ional distance ~ n9 (§2.3). 

A fundamental difference in the use of spherical coordi- 
nates is that while in a thin disk setting the vertical pressure 
gradient is balanced by the component of gravity along l z , 
gravity plays no role in balancing pressure gradients along lg 
(ALP). Gravity's place is taken by the (inertial) centrifugal 
force which for rotation on spherical shells is cx sin 2 ?9 (with 
•d = it /2 — 9). In order to compare this with "gravity along 
1 2 ", one notes that in the standard thin disk approach g z is 
oc sini? and so the fractional difference is of 0(Qn) ~ sinr?. 
It follows that the gradient in the inertial force along lg is 
gentler than that of gravity along l z (at least for $ < 1)^]. 
Note that although this somewhat diminishes the effects of 
buoyancy and the impact of the inertial force gradient (not 



5 A general relativistic analysis of this issue is complicated by the 
lack of a precise definition for the centrifugal force in a stationary metric 
(although such definition exists for a static (Schwarzschild) metric as 
demonstrated through the formalism of the 'optical reference geometry', 
Abramowicz, Carter & Lasota 1988, Abramowicz 1990, Abramowicz 
Prasanna 1990). Note, however, that the meridional hydrostatic 
milibrium equation involves an effective frequency of oscillation, Eq 
(next footnote), while the condition of radial force equilibrium, Eq 
involves deviations from Keplerian rotation rates instead (Appendix 
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gravity) on the Parker instability, for Tig /r<l this effect is 
not severe. 

Another difference is that the scale height in spherical co- 
ordinates is curvilinear. The scale height enters the general 
dispersion relation of the Balbus-Hawley (1991) instability 
by establishing a characteristic lengthscale for the (predomi- 
nant) vertical-gradient-dependent part of the Brunt- Vaisala 
frequency (c.f. their Eq [2.8]). This gradient is approxi- 
mately inversely proportional to the pressure scale height 
for adiabatic changes in the Boussinesq approximation. In 
the 'horizontal regime' k r /kg — > 0, however, this dependence 
drops out of the dispersion relation (BH92a) while simulta- 
neously yielding the maximum growth rate. The dispersion 
relation is simplified further by assuming rotation on cylin- 
drical shells (i.e. coinciding isobaric and isochoric surfaces). 
An equivalent assumption in spherical coordinates for a slim 
disk is rotation on spherical shells and with this, BH91 anal- 
ysis is straightforwardly carried on to a thicker disk. 

Non-axisymmetric Ai-type instabilities of a purely toroidal 
field only require (§[D|) that 4AS1 + k v v^ lt < for 
non-relativistic Keplerian rotation in the horizontal regime 
(Foglizzo & Tagger 1994, hereafter FT94). Although the 
geometry imposes a lower bound on k v — 1/r for nearly Ke- 
plerian rotation, the above restriction is independent of scale 
height and constraints only the toroidal field strength. To 
first order, general relativistic corrections (§[a|) do not mod- 
ify this conclusion. Indeed, magnetorotational instabilities 
in the horizontal regime are generally insensitive to vertical 
gradients and the pressure scale height enters only indirectly 
when the Alfven speed is scaled with respect to the local 
sound speed of the medium (see the discussion at the end of 



3.4) 



Lastly, since buoyant motions are commensurate on strat- 
ification and since radial pressure gradients force rotation 
on nearly spherical shells when the shear is strong (i.e. in 
a mildly advective disk), the B+ flux should naturally rise 
along the lg local pressure gradient. This is a crucial point 
that suggests that a fraction of the magnetic flux generated 
by the dynamo is promptly collimated into the funnels above 
and below the disk's surface thus creating favorable condi- 
tions for the formation of a Poynting jet. 



3.3 The dynamo equations 

The shear flow generates B v from B r on the shear timescale 
t s . In the absence of a propitious field topology for the 
development of Parker's undulate instability (e.g. FT94), 
the buoyancy of the B+ flux, which escapes the disk on a 
timescale t b , limits the growth of B v and B r and may gen- 
erate Bg albeit inefficiently: The flux buoyancy acts as an 
interchange instability whereby inhomogeneities along a field 
domain cross section induce a gradient in the buoyant veloc- 
ity along the direction of the mean field thus generating a 
small Bg component. Inclusion of this rather uncertain pro- 
cess for meridional field generation adds considerable math- 
ematical complexity to the heuristic equations and will be 



overlooked in the name of simplicity. Note, however, that 
although the meridional field is unlikely to attain strengths 
that could explicitly alter the proposed dynamo equations, 
its presence does alter the equilibrium turbulent state since 
the poloidal MRI is more efficient at feeding the turbulence 
(E. Vishniac, Priv. Comm.). This affects the IDD explicitly 
through our estimate of the magnetic Mach number for the 
turbulence (§[3~9|). 

Non-axisymmetric Af-type instabilities that develop in 
a purely toroidal field (FT94, Foglizzo & Tagger 1995 
[FT95], Terquem and Papaloizou 1996 [TP96], Papaloizou 
and Terquem 1997) contribute to both B v and B r field gen- 
eration but their growth rates decrease with field strength, 
c.f. Eq |jl7), and are likely to be important only for B r since 
shear is the primary process for azimuthal field growth. TP96 
note that because the toroidal MRI is essentially local in az- 
imuth (at least for small fields) the field need not be entirely 
toroidal, only local toroidality is required. 

Additionally, reconnection of all field components proceeds 
at a rate that is proportional to the upstream component of 
the Alfven velocity, VAif = B/-y/47rg, times the square of 
the magnetic Mach number of the turbulence (Lazarian & 
Vishniac 1999) as measured in the LRF. 

In Lagrangian, LRF coordinates, the phenomenological 
two dimensional dynamo equations are self-evident 



d t B v 
d t B r 



io _ 1 o _ 1 o 

JDy -Dip -D (p 

T s T B T r 

±B v -±B r -±B r 



(9) 



Within the heuristic phenomenology of these equations, 
when competing gain (loss) timescales differ greatly one may 
safely drop the weaker terms. 

3.4 Relativistic scalings 

In a classical Keplerian disk, the timescales related to shear, 
magnetorotational and local vertical oscillations of test par- 
ticles are all scaled by the local Keplerian rotation frequency 
(as meassured by observers at rest in a global frame). In a 
relativistic slim disk setting, however, the timescales related 
to all of these processes are slightly different. The importance 
of this remark becomes evident when one attempts to gen- 
eralize the instability scalings (LRF time and length scales) 
for a slim disk in the Kerr-Newman geometry. 

Through Oort's A parameter (see below), the shear and 
magnetorotational timescales are intrinsically related to the 
azimuthal rotation frequency of the disk. For nearly Keple- 
rian prograde rotation the natural scaling is simply 51+ (c.f. 
Eqs @ & Q). 

On the other hand, the natural frequency of local merid- 
ional oscillations of the four velocity GP98^|, Cl z , which is di- 
rectly involved in determining the local pressure scale height 

6 obtained by expanding Carter's (1968) fourth constant of geodesic 
motion in the Kerr geometry, Q, about i? = it/2 — 8 and taking an affinc 
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Tig, c.f. Eq ||, differs from the prograde Keplerian rotation 
rate by a factor 5ft, where (Riffert & Herold 1995) 



5ft 2 = ( ^ 



1 - 4ar~ 3 / 2 + 3aV 



(1 + ar- 3 / 2 ) 2 (l - 3r-! + 2a7- 3 / 2 ) ' 



(10) 



As advertised, the growth rates of Al-type instabilities 
are characterized by a relativistic generalization of Oort's 
A parameter while for buoyant motions and reconncc- 
tion timescales (and P-type instabilities) the characteristic 
timescale is (loosely speaking) the meridional 'Alfven transit 
time': Tte/v^g. Thus, although the pressure scale height, 
Tis, naturally normalizes all meridional lengths, similarly 
scaled magnetorotational wavenumbers must be corrected by 
a factor 5R from their classical counterparts 



• Rel 



Sftfc 



Class 



(11) 



That this is indeed the case may be more easily seen by 
normalizing either the classical (with fi = i v /x r ) 



4A ort^+(k- VAlf) 2 < 



3) 



or the generalized relativistic stability criterium for shearing 
instabilities, Eq p6| ]. Because for prograde Keplerian orbits 
the term driving the radial differential force scales as 17+ , 
normalization of the Alfven velocity to the local sound speed 
in the (k • vam) 2 term yields the correction factor if one de- 
sires to normalize the wavenumber to the disk's pressure scale 
height. 

3.5 Linear shear 

We choose to model the shear flow in a linear shearing sheet 
approximation (Goldreich & Lynden-Bell 1965) noting the 
well known limitations of such an approach, i.e. ignoring 
global curvature (Ogilvie and Pringle 1996, TP96) and non- 
linear shear effects. The influence of shear is thus parameter- 
ized by a relativistic generalization of Oort's first constant, 
t^ 1 = 2Ag' a I s t s = d\ nr Vt (= | £!kc P , for a Keplerian disk). 

A relativistic expression for Oort's A value that contains 
the radial gradient of I and (3 r implicitly is (Appendix §|a]) 



af<p(r, d r n) = Aj 2 



1 d In r 



A Rcl 
A Oort • 



This result differs subtly from the one in Novikov and Thornc 
(1974) where 7 goes to "f v for a thin disk with negligible radial 
flow. 

For nearly prograde/retrograde Keplerian orbits, from Eq 
l35| one finds 





1 










derivative 










Q = Ul 4 



n± 



1 ± ar- 3 / 2 



9 a Rcl 



a 2 (S 2 - 1) 



(12) 



Shear-fed instabilities, on the other hand, pertain to growth 
rates slightly different from due to relativistic effects^. 

Shear also forces the radial wavenumber of perturbations 
to evolve according to 



k r (t) 



2 A_/c^3 1 . 



(13) 



However, we shall concern ourselves only with the horizontal 
regime of magnetorotational instabilities, k r /kg — > 0, noting 
that this constraint yields the highest growth rate from both 
the Balbus-Hawley (BH92a) instability (§|3.6|) and toroidal, 



non-axisymmetric A^-modes (§3.7) 



3.6 The poloidal (Balbus-Hawley) MRI 

Classically, the axisymmetric magnetorotational instability 
(Balbus & Hawley 1991) only requires that d r ft < and 
that the Bg field be "weak" in the sense v^ li < (1/tt) \Aft\ H 
which is derived by imposing an upper limit (27i) on the 
'vertical' wavelength in the Keplerian thin disk limit (note 
that the radial field component is irrelevant in the horizontal 
regime: k r /kg — * BH92a). 

On the other hand, the classical instability criterium 
(Chandrasekhar 1961 and Appendix §0), 4Af2 + fc e u e Alf < 0, 
sets a upper limit on kg 1 ** = \/3D,/v^ u beyond which the 
energy cost of bending field lines on small scales suppresses 



,Aif 



the instability. The fastest growth occurs for kg = 
a rate that equals Aq 1 ^ as long as kg ^tt/Hq (- 

n/ Cs ) 

Although a fully general relativistic stability analysis of 
MHD disk systems in the Kerr geometry is beyond the scope 
of this paper, in Appendix §^|wc derive an approximate (see 
footnote 6) relativistic timescale, corrected for the form of the 
"radial differential force" in the Kerr geometry. For nearly 
Keplerian orbits, this procedure leads to replacement of the 
term driving the (classical) radial differential force d\ nr il 2 



by 4f2-|-A^ ± For prograde orbits: 
estimate the timescale and Eq [|lT| 
wavenumber, one finds 



using Eqs |3J & |37j] to 
to normalize the optimal 



1 




3 2 


'A' 


T R.cl 

BH 


max 


4 ' 


V_ 



(=-A 



eff 



(compare this with Eq JT^]), and 



■y 1, - /2 c s ^ Sl+ _ 1 2 c s 



r v, 



(14) 



(15) 



We identified the r.h.s of Eq p4[ with the maximum growth 
rate of shear-fed instabilities (such as A'l-type) in relativistic 
prograde Keplerian Kerr disks as meassured in the LRF of 
the fluid: 

I = _ A eff 



Rcl 



max,+ 



7 The relativistic shearing and magnetorotational timescales are de- 
rived by modifying only the relativistic analog to the classical "radial 
differential force" (FT94), and assuming that, both, the centrifugal and 
the coriolis forces act as they do classically, c.f. Eqs ]3p| and footnote 
and discussion at the end of Appendix aAl 
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3.7 The non-axisymmetric toroidal MRI 



These estimates are in fact subject to the condition 



Although BH92b originally found that the most unstable 
magnetorotational modes correspond to axisymmetric modes 
evolving in a poloidal field configuration (through desta- 
bilization of Alfven modes), several authors (FT94, TP96, 
Matsumoto & Tajima 1995) find similar growth rates for 
non-axisymmetric modes in a purely toroidal field configu- 
ration (albeit at larger poloidal wavenumbers, Balbus 1998). 
Note that in hyper-accreting 7-ray burst black holes, non- 
axisymmetric modes may be strongly excited if these feed 
the angular momentum loss through an associated gravita- 
tional wave emission. 

This instability corresponds to destabilization of the slow 
MHD mode of wave propagation (FT94, TP96) as is also the 
case for Parker's undulate instability. FT94 (see also TP96, 
Shu 1974) remark that the inertial effect of rotation enters 
through the coriolis force term, slightly twisting flux tubes 
of displaced fluid elements in the azimuthal direction, and 
through the radial differential force. 

If the ratio k r /ke <C 1 (=>• 'horizontal' regime), meridion- 
ally localized disturbances are more prone to the radial dif- 
ferential force. If this force is strong enough to overcome the 
magnetic tension this results in the non-axisymmetric mag- 
netorotational instability (FT95). Radially localized pertur- 
bations kg/k r <C 1 (=> 'asymptotic' limit, c.f. Eq 13|), on 
the other hand, would naturally evolve into Parker unstable 
modes if the field topology were favorable. More generally, 
this instability generates both B r and B v (in the horizontal 
regime, the slow mode is essentially an Alfven wave) and af- 
fects the field configuration by twisting of toroidal field lines 
into flux ropes perhaps enhancing turbulent pu mpin g (Vish- 
niac 1995a) and the associated field buoyancy (§ |3.8|) . 

We use the Kerr geometry generalization of the radial dif- 
ferential force, §[a| and define a v = PB^/(Pgas + Ph) and a 
rescaled tm and A°5 (< 0) 



- A < 1 + (2a + 1) 



(19) 



A eft 



(16) 



to adopt the results in FT95 as approximate relativistic gen- 
eralizations of the toroidal MRI scales. 

For prograde, nearly Keplerian rotation 

, A cff ,2 2 _ 



(1 + 2a»)(l + 2aJl + A)) + 1 + a v (2 + A) \ (17) 



and 



2 / c. 



r y ,,\„ l >< (2A+(l + a v )f M 2 ) 



x S 



r v v Ali 



(18) 



In a stronger effective shear flow the optimal wavenumber in 
Eq |l^] becomes imaginary and the toroidal field is unstable 
to a faster radial interchange (k v = 0) instability according 
to the Rayleigh criterion (T. Foglizzo, Priv. Comm.). Inter- 
estingly, this destabilization process is mathematically simi- 
lar to the one that occurs when suprathermal particles push 
the Parker instability into interchange-like mode (FT94) for 

Psuprathcrmal/Pgas > 3 + 4a (see below). 

To the IDD, this means that when condition |l9| is not 
meet, the toroidal MRI does not grow the radial field nor 
does it provide structure nor field reversals in the linear anal- 
ysis. Moreover, note that when the wavenumber becomes 
very small, the local analysis is out of its intended regime and 
global curvature effects become important. The IDD is thus 
inadequate in this limit. Fortunately, condition p9| is only 
meet at most within r§5 [GM/c 2 ] for a non-rotating black 
hole and it moves inward for a rotating black hole (§4.2). 



for the non-axisymmetric MRI modes of most vigorous 
growth. 



3.8 On the buoyancy induced by 
radiation and e pairs 

One may physically expect that pairs would enhance the 
buoyancy of magnetic flux if such a lighter fluid component 
coupled to the 'other' light relativistic fluid of the medium, 
i.e. the magnetic field. Yet, further elaboration on such 
a conjecture requires some understanding of the dynamics 
of quasi-equilibrium pair plasmas in magnetic fields well in 
excess of the Q.E.D. field scale (Bq = m 2 /e(c 3 /h) ~ 44. 
TG), a problem that has yet to be explored. 

For instance, under this circumstance pairs are created 
(annihilated) mainly via I7 decay(production) with strongly 
asymmetrical energy profiles for each member of the pair 
(Daugherthy & Harding 1983, Harding 1986). It is thus 
conceivable that at very high disk temperatures the pairs 
may have considerable energy density in a non-thermal dis- 
tribution of suprathermal particles. This being the case, 
Parker's suprathermal mode of hydromagnetic wave prop- 
agation (which strictly speaking applies in the collisionless 
regime) will enhance Parker's undulate instability possibly 
pushing the latter into an exchange-like mode (FT94) with a 
normalized growth rate (wcrowth/^) 2 oc (a± — /3)/(l + 2a±), 
where a± is calculated from the Bj_ field, i.e. including B r 
in contrast to Eqs |l7| k. |§. 

On the other hand, (thermal) transrelativistic pairs are 
well coupled to the nucleon component even when coulomb 
dragg is low (due to the neutronization of the material) 
and ambipolar diffusion of the field+pair component through 
neutron matter is dismal (although the turbulence may help 
the latter). Thus, the most likely mechanism for fast horizon- 
tal flux escape would have to involve nucleon density deficits 
inside flux ropes or, in the case of a diffuse field, turbulent 
diffusion. 
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A plausible mechanism that promotes baryon unloading 
from field lines is turbulent pumping (Vishniac 1995a) by the 
MRI which must favorably compete with turbulent diffusion 
of matter back onto flux ropes. Under this assumption, the 
stretch, twist and fold of field lines by (enthalpy-weighted) 
sub-Alfvenic turbulence augments the field energy density 
and releases matter from field lines that would otherwise be 
"frozen-in" . For the marginal case of Alfvenic turbulence^, 
nearly empty B± flux ropes in a gas pressure dominated 
disk acquire a drag limited buoyant velocity vt cx {v^) 2 jc s . 
Moreover, assuming efficient diffusion of radiation and 
pairs into the flux tubes, in this picture the buoyancy loss 
rate, Vb/Tte, is enhanced by a factor § p/p ga , s (Vishniac 
1995b). 

On the other hand, in a diffuse field configuration the local 
value of the turbulent helicity implies migration of the field 
to flux poor regions of the accretion disk (E. Vishniac, Priv. 
Comm.). Phenomenologically, this is also born in numeri- 
cal simulations with vertical stratification (Brandenburg, et 
al. 1995, Stone, et al. 1996). Since the main contribution to 
the field generation occurs in regions of highest pressure, i.e. 
the disk midplane, this argument yields a systematic merid- 
ional motion of the field (up to a gradient in the diffusion 
coefficient) 



2 ^tnrh^cor 

Vh = 

r He 



r c s 



mi 



(20) 



where r cor ~ is the correlation timescale of the turbulence 
and Mg = w turb / w Alf [ s its magnetic Mach number on the 
largest eddy scale. 

The estimate for the field loss rate from the disk depends 
on whether or not the radiation/pairs enhance the buoy- 
ancy loss rate. In the context of a diffuse field configuration 
(which is on more solid grounds), this enhancement is only 
a conjecture. With this in mind, using Eq |2(| and defining 
£ = Pg as /p, we write the characteristic time as a meridional 
Alfven 'transit time' (admittedly a misleading nomenclature 
given that Alfven waves move along field lines) 



K 

Vb 



(B v , B r ) = L2L= v *= r f[J2 s \ 

\ Si a I "4- 



c. 



where 




(21) 



(22) 



3.9 Turbulent reconnection 

It is well known that in spite of generally predicted low re- 
connection rates (Sweet 1958, Parker 1979), well studied sys- 
tems such as the solar corona and chromosphere indicate 
that reconnection is fast once it begins (Dere 1996, Inncs 

8 Note, however, that numerical experiments bias our expectation of 
MRI turbulence toward super-Alfvenic values, i.e. the Maxwell stress 
tends to dominate over the Reynolds stress. 



et al. 1997), essentially occurring at an order of magnitude 
below the Alfven speed. 

Following this observation, Lazarian & Vishniac (1999) 
show that in the Goldreich & Sridhar (1997) model of strong 
MHD turbulence, reconnection in a weakly stochastic field 
occurs at a fraction of the Alfven speed which equals the 
square of the magnetic Mach number of the turbulence M^ . 
In a medium where MRFs feed the turbulence, an educated 
guess of this number is (Vishniac & Diamond 1992, Zhang, 
Diamond & Vishniac 1994) 

Mfe[(/ lf / Cs ) 1/3 , l.]. 

We adopt the lower value noting that the Balbus Hawley in- 
stability is likely to be more efficient at feeding the turbulence 
and that numerical simulations generally show a dominant 
Maxwell stress regardless of the initial seed field. 

In keeping with our informal nomenclature, the reconnec- 
tion rates are written as inverse Alfven transit times calcu- 
lated from the component of the field that undergoes recon- 
nection 

^ L rcc _ _ _ rec Q — 1 

l^ l c is a perpendicular lengthscale associated with the mean 
distance for reversal of the field component undergoing re- 
connection. Let us label according to 



Sx r 



X, Sa 



V, and 5x° 



Z. 



with a subscript denoting the component of the field under- 
going reconnection. 

The perpendicular lengthscales associated with the mean 
distance for field reversal are derived from the fundamental 
linear lengthscales for coherent field pumping. The latter 
are supplied by half of the toroidal MRFs (wave)length scale 
(recall that in the horizontal regime one deals essentially with 
Alfven wave destabilization). 

Azimuthal, radial field reversal, Y r , directly involve s th e 
optimal wavenumber of the non-axisymmetric MRI (§3.7). 
This lengthscale is sheared in the azimuthal direction thus 
coupling it to the radial, azimuthal field reversal lengthscale, 
X v . Note that the meridional lengthscales associate with 
the poloidal MRI are not relevant since the meridional field 
is assumed to be weak (albeit quickly unstable). 

It has been argued (TP92) that the shear process generally 
reduces radial lengthscales until these supply the fundamen- 
tal reconnection channel. This is indeed the case when any 
of the coherent field pumping process is slow when compared 
to shear, r 3> t s ~ fiT ; but this is never the case when the 
MRI is the culprit (even for the slower manifestations of the 
instability). One may thus considerably simplify the calcu- 
lation by ignoring this aspect of shear in the reconnection 
process. 

On the other hand, following TP92, we obtain an esti- 
mate of the radial, azimuthal field reversal lengthscale, X v 
by noting that the time evolution of wavenumbers implied 
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by Eq jll| during one MRI timescale couples the azimuthal 
lengthscale to the radial lengthscale, i.e. ly — (t m /t s ) x . 

Thus 



7T / T \ 7T 

Y r = — and X v = f f- 1 x — . 



(23) 



3.10 Equilibrium Solutions 



In general, Eqs |9j] are implicitly complex due to the 
non-linearities introduced by the field dependence of gain 
timescales, Eqs |l7j, the buoyancy loss rate Eq 
and the reconnection times which involve the perpendicular 
lengthscales Eqs [ p3| . 

Scaling wavenumbers to the inverse pressure scale height 
k — > k/H.Q, and defining 



and 



He 



A 



.4^ = f- and W- £ ' = ^ = 41f- (24) 



A 



We 



A' fc. 



and 2A' 



^5. 

n 4 



(c.f. Eq |l6|], although these are positive definite), a set of 
normalized dynamo equations for prograde orbits follows 



dt>B' 



2A'B' r -^B' ip B' ± -^ r (B' ip r, (25) 



1 



1 



d v B' r = AB' lp --B' r B' ± --(B: 



>r 



(26) 



where the fields are in velocity units and normalized to the 
soundspeed (£>' = B x y/4:TTgc s ) and the time is normalized 
to the inverse of the prograde Keplerian frequency t' = fl + t. 

In a steady state, these equations must satisfy 



BL 





[Mf ] 2 x < 


f 5R 2 


2A' 




K 


[Mf r x | 




A 





7T A J 

7T BL 



(27) 
(28) 



which we solve using a multidimensional Newton-Raphson 
method (e.g., Numerical Recipes in C, Press et al. 1988). 



4 A toy calculation 
4.1 Self-critique 

The instability driven dynamo depends self-consistently on 
the relativistic shear, 2A', and MRI, A, timescales, on the 
coherent field pumping lengthscale supplied by the MRI, 5, 
and on the magnetic Mach number of the turbulence, Mg. 
A phenomenological scaling for the buoyant velocity in the 
turbulent medium, Eq 20 1, follows from the limiting cases of 
nearly empty flux tubes and turbulent diffusion of a spread 



out field (and by analogy to other buoyant instabilities such 
as the Parker instability). 

An open issue is whether a radiation pressure dominated 
environment helps the buoyancy of the field and what this 
means for the magnetic energy deposition rate, particularly 
in the context of a 7-ray burst. If the emerging magnetic 
flux (in ropes) is relatively baryon-free, this helps its escape 
from the disk. On the other hand, the meridional diffusion 
of B^ flux in a semi-thick disk setting also tends to separate 
the field from the (heavier) baryon component, aiding in the 
escape of flux -(-radiation and pairs. 

In order to address the overall energetics and to gain a 
qualitative feel of the postulated hydromagnetic energy con- 
version process, we have constructed a toy model which 
adopts the pressure ratios from the standard model of PWF 
and (non-self-consistently) assesses the energy loss from the 
buoyancy of the field. Evidently, a realistic model should 
account for the back-reaction of the field on the flow which 
in the context of an accretion disk means primarily the an- 
gular transport implied by the -rip- component of the mag- 
netic stress (note that in the notation of the previous section, 
a™ s ag = B' r B' v ). In the absence of other significant sources of 
'anomalous' viscosity, a self-consistent hydromagnetic accre- 
tion disk model is possible in principle when our equations 
are solved with a compatible set of relativistic hydrodynamic 
equations (e.g. GP98, Abramowicz et al. 1997) under the as- 
sumption that the angular momentum is transported by the 
largest magnetic eddies. 



The radial flow (§3T) is passively accounted for by carry- 
ing out the analysis in the comoving frame and assuming that 
the linear timescale of the non-axisymmctric MRI faithfully 
reflects the growth rate of the radial field (an assumption that 
is not free of criticism). On the less optimistic side, the radial 
interchange instability (§[3jj) of the Bj^ flux in a strong effec- 
tive shear field shuts the dynamo off at the innermost section 
of the disk where most energy could be extracted. This in- 
stability is already present when 1 < —A < l + (2a + l) _1 (T. 
Foglizzo, Priv. Comm.) but with a slower growth rate than 
the MRI. Thus, even in the absence of substantial meridional 
flux (Park & Vishniac 1994), radial buoyancy may change 
the picture we have portrayed hereby (on the brighter side, 
the radial interchange opposes the advection of the magnetic 
field). In this regard, note that condition |l^| depends only 
on the toroidal field strength when k v — > but we cap the 
latter at k v = 1/zu 



4.2 The magnetic energy deposition rate 

The local magnetic energy flux (from buoyancy) is obtained 
from the part of the time derivative of the energy density 
which is lost to buoyant motions 



Ti-efs — «b = 



4vr I 



GM 



-3/2 I> £ [M B ]2 x Bl 

£ A 



(29) 
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where B r and B v equal the second terms in Eqs [£5] & |2FJ . 

Ado pting the results from the standard model of PWF (c.f. 
§§ & 3.1), this form for the specific power output is used 



in four estimates of the steady magnetic flux deposition rate 
for r £ [4.5, 40]: non-rotating black hole with (hBuoy-a.00), 
and without enhanced buoyancy (lBuoy-a.00), and nearly 
maximally rotating black hole, a = .95, with (hBuoy-a.95), 
and without (lBuoy-a.95) enhanced buoyancy. Curiously, be- 
cause the fields are higher in the low buoyancy scenario, these 
yield the highest output rates (in erg sec -1 ): 8.I+51 (lBuoy- 
a.00) and 1.0 +52 (lBuoy-a.95) vs 5.1 +5 i (hBuoy-a.00) and 
6.6+51 (hBuoy-a.95). 

The 'half-luminosity' radius displays some interesting be- 
havior as well: for hBuoy-a.00, the IDD becomes opera- 
tional (c.f. radial interchange instability) at r m i n ~ 5.3 and 
rc — 10.3, while for lBuoy-a.00, r m ; n ~ 5.9 and rc — 9.2, 
because of the higher stationary fields in the low buoyancy 
case: [B' v , B' r ] = [.72, .80] vs [.42, .48] at r min , respectively. 

For nearly maximally rotating holes, the IDD operates 
down to smaller radii because frame dragging yields a slower 
effective shear: r m i n ~ 5.4, B' = .93 and rc — 8.5 for lBuoy- 



a.95 and 7',, 



4.6, B' = .57 and r c ^ 8.5 for hBuoy-a.95. 



The implied ass-p^ameter at the minimum operational 
radius is highest for maximally rotating, low buoyancy disks 
« .80, and smallest for non-rotating, highly buoyant disks 
w .21. This magnetic stress decreases quickly at larger radii: 
for low buoyancy, ass ~ -39, and .31 at r = 10, and 20, 
respectively; while for high buoyancy, ass ~ -16, and .17 at 
r = 10, and 20. 



fects, the field generated by the n on-a xisymmetric MRI in 
a strong effective shear flow, c.f. §3/7, pushes the optimal 
wavenumber to its lower bound which implies pumping of 
coherent field structures with a lengthscale of 0(ttuj) ~ nr. 
In other words, the coherence lengthscale of the toroidal MRI 
in a strongly sheared rotor is large. Thus, in principle all the 
energy deposited at the innermost radii may go into a Poynt- 
ing flux. In the low buoyancy, high spin case, our estimate 
indicates that this may be as much as 10% of the total output 
rate, without invoking inverse cascade arguments. Invoking 
the Blandford-Znajek mechanism can only increase the above 
estimate. If fully self-consistent models yield similar results, 
this could render the proposed hydromagnetic energy conver- 
sion mechanism a more viable alternative to neutrino-burst 
driven models of 7-ray bursts. 
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4.3 Poynting jets and 7-ray burst engines 

We have constructed a reasonable set of heuristic two di- 
mensional dynamo equations for magnetic field components 
in the comoving frame under the premise of negligible gen- 
eration of meridional field. These equations approximately 
account for field generation by the shear flow and by the 
non-axisymmetric MRI in the Kerr-Newman geometry. Two 
local processes are invoked as field loss terms: turbulent flux 
buoyancy (or vertical diffusion) and turbulent reconnection. 
Self-consistent equilibrium solutions to these equations are 
found in proximity to equipartition field values by adopting 
pressure ratios and the scale height from the standard model 
of PWF. Although, one must not take these estimates at face 
value until a self-consistent calculation is displayed, the to- 
tal hydromagnetic energy deposition rate from buoyancy is 
found to be at least comparable to the neutrino luminosity 
in that scenario, L v ~ 3. 3+51 erg sec -1 . 

The magnetic field is removed from the disk interior and 
its ultimate fate depends on the details of the disk corona 
dynamics which I do not address here. At 9 ~ ®h-i some 
of this flux can be expected to undergo an inverse cascade 
to form larger coherent field structures (Tout and Pringle 
1995). This large scale field is the relevant one with regards 
to Poynting jet production. However, ignoring curvature ef- 



A Relativistic Shearing Instabilities 

In a local kinematic analysis of the Hill (1878) equations 
for a "sticky disk", Balbus & Hawley (1992a) have shown 
that the equations of motion for a (corotating) Lagrangian 
displacement vector £ (with l r — > l x , and l v — > l y ) 

C !S 20$ = -{4A 00 rtfi + X: x K a 

C« + 20q = -Ky i". (30) 

yield a dispersion relation identical to the one derived from 
an MHD stability analysis of axisymmetric perturbations 
of a purely poloidal disk field configuration in the limit 
k r /k0 — ► 0, provided that 

K-X = fcy = kgV&tfg. 

Notably, although both terms in the curly brackets of Eqs 
are stabilizing, their combined influence is destabilizing 
(Balbus & Hawley 1996) because of the spatial coupling. 

Foglizzo & Tagger (1995) find an identical set of equa- 
tions for the magnetorotational instability of a purely toroidal 
field in the horizontal regime, provided that the "spring con- 



11 



stants" are replaced by 



and /C„ 



2 2 
C s v Alf,ip 



C 2 + V 2 



AK,<p 



This is not surprising given the generality of the magnetoro- 
tational destabilization process (BH92a). 

In either case, the dispersion relation is of the form 
t- 4 - br- 2 + c = 0, with b = 4(0 2 + AO) + QC X + JC y ), and 
c = /C„(4Af2 + IC X ), where we have set d\ nr fl 2 = 4AS1 (recall 
A-Oorif = hd]nr&)- This dispersion relation yields a general- 
ized classical instability criterium (i.e. t 2 < 0) 



4Ao ort n + (k • v A i f f < 



(31) 



for magnetorotational instabilities in either poloidal or 
toroidal magnetic field topologies, in the horizontal regime 
k r /ke — > 0. Because in most cases of interest to us Aoort < 0, 
the first term in Eq |n| drives the destabilization process. 

Classically, this term (4Af2 = d mr fi 2 ) may be derived by 
considering the unbalance of inertial forces induced by de- 
viations from Keplerian rotation of radially displaced fluid 
elements, i.e. induced by the imbalance of centrifugal vs. 
radial gravity (FT95). Alternatively, in the Hill equations 
(which follow from Galilean transformations of the Newto- 
nian equations of motion) the equivalent expression, which 
contains Oort's A constant implicitly, may be obtained by 
expanding the two factors in — "r f2 2 "- about the equilibrium 
point and keeping only linear terms in Sr. 

In the general theory of relativity, on the other hand, the 
effect of gravity is geometrically included in the equation 
of motion and the meaning of force is intrinsically different 
from that in the Newtonian theory even if the acceleration is 
meassured in a frame instantaneously at rest (Abramowicz 
& Prasanna 1990, Abramowicz 1990). Still, in the LRF a 
procedure similar to the radial differential force analysis (c.f. 
Eq [35-37] of FT95) may be used to derive an approximate^ 
relativistic analog to Eq [|TJ. 

Although geodesic Keplerian flow does not 'feel' gravity, 
deviations from relativistic Keplerian rotation unsettle a flow 
in otherwise radial equilibrium (ignoring radial pressure gra- 
dients). The equation that describes this behavior was first 
derived by Lasota (1994, see also Abramowicz et al. 1997 and 
GP98) 



ri- 



(32) 



Recall that (3 r is the radial speed as meassured by an observer 
at rest in the corotating (CRF) frame. 

9 Eqs fcifj] are derived under the premise that the standard centrifugal 
and coriolis force terms act in the classical sense. In the Kerr geometry, 
a precise definition of the centrifugal force, as formally required to de- 
rive a relativistic analog to Eq Jjl) , is hampered by the appearance of a 
coriolis-type radial force related to the Lenz-Thirring effect (Abramow- 
icz, Carter & Lasota 1988). These forces, which are not included in our 
approximate analysis, are important at the innermost region, 3SrS5, 
of accretion disks in the Kerr geometry (Chakrabarti & Prasanna 1990, 
Prasanna & Chakrabarti 1990). Their combined effect will be explored 
elsewhere. 



Noting that the radial speed as meassured in the LRF by 
definition equals zero (albeit instantaneously) we re- write the 
l.h.s. of Eq (32) in terms of the radial gradient of this v r 



d r 0l = — 7 r 4 d r v 2 . 



(33) 



Since the gravitational potential is implicit in the geom- 
etry, the r.h.s. of Eq |32| may be interpreted as a New- 
tonian acceleration meassured in the local rest frame, aka 
p r = —d r H = —d r T for a "free" particle (note that the 
LRF is the only frame where, loosely speaking, it is safe to 
treat a 'force' as a having its Newtonian meaning). 

Next, we expand the two terms in parenthesis using the 
definitions for the prograde and retrograde Keplerian fre- 
quencies in BLF coordinates ft± = ±(r 3 / 2 ±a) _1 and Taylor 
expand SI 2 around either Keplerian frequency 



1 



n 
si. 



1 



{1 - 2aVt 



n 2 (r 3 - a 2 )} 



a 



where the difference Sfl 2 is taken with respect to (frame- 
dragged) relativistic Keplerian rotation rates. 

Lastly, putting together Eqs g3| |||] to construct the 
radial differential force, defining the relativistic shear rate in 
the LRF 

Vll^z) (35) 



aRcI 
A Oort 



d\ni 



(recall this expression contains the radial gradients of £ and 
[3 r implicitly); and noting that all other quantities in Eq 
pit are locally defined values, the desired relativitic stability 
criterium obtains 



" 1 " 




V 


r 3 ( 



1 



a. 



T7r - Ag c o 1 rt O ± + (k.v Alf ) 2 <0. 



(36) 

(recall 0_ < 0). 

The timescale derived from the relativistic modification to 
the radial differential force follows straightforwardly by re- 
placing the expression in the curly bracket of Eqs j3(]] by its 
relativistic equivalent in the l.h.s. of Incq fjil , and following 
the analysis of BH92a thereon. In the limit k r /kg — ► (the 
horizontal regime) this timescale should be identified with 
the maximum growth rate for (magnetic) shearing instabili- 
ties 



1 




-7 2 


'A' 


A dVL \ 


T Rcl 
M 


max 


V 


\ 2 d\nr J 



l{r 3 -a 2 + a (ri±a)}. 



(37) 

This expression reduces to the correct limit when a — > 0, with 
the shear force reversal at r — 3 embodied by the expression 
for Ao^t as remarked by GP98 (see also, Anderson & Lemos 
1988, Abramowicz & Prasanna 1990). 

For Keplerian, prograde orbits; using Eqs (TT 
finds 



one 



A' 



off _ 



1 


3 2 


'A 


T Rcl 

M 


= 4 7 

max 


V 



I!. 



(38) 
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while for retrograde disks, the strength of the instability is 
somewhat diminished from its classical growth rate. 

In the horizontal regime, one may expect the maximum 
growth rate of shear-fed instabilities (above) to differ from 
the relativistic shear rate (compare Eqs [|l4| & [l^]) because 
the former corresponds (classically) to obliquous local fluid 
excursions with £ x = —£ y (see §2.4 of BH92a) while the 
gradient of the velocity field (i.e. the shear) is maximized 
for excursions perpendicular to l y . Since the local thetrad 



of basis LRF vectors is rotated by the geometry (c.f. §2.1), 
it is not supprising that new metric corrections are involved 
in the growth rate for magnetorotational instabilities. While 
this argument does not prove the validity of Eq |37j (which 
is, after all, an approximate result), it does motivate our 
finding. 
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